Decoding mode-mixing in black-hole merger ringdown 
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Optimal extraction of information from gravitational-wave observations of binary black-hole coa- 
lescences requires detailed knowledge of the waveforms. Current approaches for representing wave- 
form information are based on spin-weighted spherical harmonic decomposition. Higher-order har- 
monic modes carrying a few percent of the total power output near merger can supply information 
critical to determining intrinsic and extrinsic parameters of the binary. One obstacle to constructing 
a full multi-mode template of merger waveforms is the apparently complicated behavior of some of 
these modes; instead of settling down to a simple quasinormal frequency with decaying amplitude, 
some \m\ 7^ £ modes show periodic bumps characteristic of mode-mixing. We analyze the strongest 
of these modes - the anomalous (3, 2) harmonic mode - measured in a set of binary black-hole 
merger waveform simulations, and show that to leading order, they are due to a mismatch be- 
tween the spherical harmonic basis used for extraction in 3D numerical relativity simulations, and 
the spheroidal harmonics adapted to the perturbation theory of Kerr black holes. Other causes of 
mode-mixing arising from gauge ambiguities and physical properties of the quasinormal ringdown 
modes are also considered and found to be small for the waveforms studied here. 
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I. INTRODUCTION 

Since the first successful simulation of black-hole bina- 
ries (BHBs) through late inspiral, merger, and ringdown 
in 2005 [l|-y| , theoretical interest has centered on the re- 
sulting gravitational waveforms. A crucial tool in wave- 
form studies has been the analysis of the radiation wave 
pattern in spherical harmonic components. This decom- 
position is useful both in the physical interpretation of 
the radiation, and in structuring the waveform informa- 
tion content for the development of approximate analytic 
or empirical encodings. 

The self-consistency of results for the dominant 
quadrupolc waveforms across numerical codes was 
quickly established [3, enabling rapid study of the 
basic characteristics of mergers p-[l3j Researchers soon 
began to build analytic template models compatible 
with these numerical results as well as with the post- 
Newtonian (PN) at earlier times, to provide relatively 
quick waveforms for specified BHB source masses and 
spins While expected to be sufficient for detec- 

tion of BHB mergers, quadrupole-only templates will not 
lock down most of the intrinsic (masses, spin magnitudes 
& directions) and extrinsic (sky position, phase) BHB 
system parameters. To gain an understanding of these 
parameters requires a richer template bank, one that in- 
cludes all of the relevant angular modes of the signal 

MM. 



Working with a spherical harmonic basis of spin- weight 
s = — 2 [2l], [2^], several studies [23H r 27j have found that 
after the dominant quadrupole [l = 2, m = ±2) modes, 
the next most important modes tend to be the higher 
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FIG. 1: Gravitational-wave luminosity from the merger and 
ringdown of a nonspinning black-hole binary of mass ratio 4:1, 
decomposed by harmonic mode. 



to = ±1 modes: (3, ±3), (4, ±4), etc., though odd-m 
modes are sometimes suppressed by symmetry. We have 
also seen, however, that certain to < \£\ modes can be 
important. Prominent amongst these are the (2, ±1) and 
(3, ±2) modes. Figure Q] shows the radiative power for 
the most important modes in the case of the merger of a 
4:1 nonspinning BHB. Here we see that the (2, 1) mode 
has actually overtaken the (5, 5) mode in importance by 
merger time. 
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A key feature of BHB mergers exposed through the 
spherical harmonic decomposition waveform studies is 
the rather clean separation of the sometimes compli- 
cated mix of signal frequencies, achieved by angular 
mode decomposition. Even when typical observers would 
measure complicated wave shapes combining several fre- 
quency harmonics, these harmonics largely reduce to 
slowly evolving sinusoids in each spherical harmonic com- 
ponent mode. To a very good approximation, this struc- 
ture holds consistently through the inspiral, merger, and 
ringdown p3l I281430T ] . This pattern of frequency separa- 
tion is extremely convenient in allowing relatively simple 
encodings of the waveform information in analytic mod- 
els. 

Partly because of these properties, angular-mode de- 
composition has become a standard approach to compar- 
ing waveform simulations with each other, with analytic 
post-Newtonian calculations, and with developing empir- 
ical waveform template models. These uses of the decom- 
position technique have elevated its significance from its 
beginning as an interpretive convenience to its current 
status as an essential component of how we quantita- 
tively understand gravitational-wave signals. Thus we 
must be aware of the possibility that artifacts of arbi- 
trary choices in the details of the decomposition proce- 
dure may interfere with our quantitative understanding 
of the waveforms themselves. 

Such concerns are particularly notable when we see 
unusual features in the decomposed waveforms seeming 
to violate the a posteriori expectation of clean separa- 
tion of frequencies. Several authors @, l25T - [27l l3l| have 
noted that the (3, 2) mode in particular typically seems 
to break from this simple pattern, showing unusual post- 
merger features that require investigation and resolution 
before a useful model can be developed. In some of the 
earliest merger simulations, Buonanno et al. Q already 
noted the presence in the post-merger "ringdown" (3, 2) 
mode of both (3,2) and (2,2) quasinormal-mode (QNM) 
frequencies. 

In this paper, we investigate these (3, 2)-mode anoma- 
lies, with a survey of 3D numerical simulations of the 
merger of various comparable-mass BHBs with non- 
precessing spins, exploring a range of possible "causes" . 
We find that the dominant part of the measured mode- 
mixing that underlies the anomalous effect can be at- 
tributed to our use of spherical harmonics rather than 
the spheroidal harmonics expected by Teukolsky pertur- 
bation theory 

The remainder of this paper is laid out as follows: In 
Section [TU we review the numerical evidence for mode- 
mixing in existing (3, 2) evolutions, and show how well 
it is captured by a simple two-mode phcnomenologi- 
cal model for the ringdown waveform segment. In Sec- 
tion [nil we discuss general models for why mode-mixing 
should be expected, including effects of coordinate dis- 
tortions in the radiation extraction spheres, and of ill- 
adapted harmonic basis functions in the radiation decom- 
position. In Scction lPVl we introduce our set of expanded 



numerical evolutions, arranged into "equivalence classes" 
of common end-state Kerr spins, which we analyze in Sec- 
tion [Vj fitting the measured contributions of two-mode 
models to our models. We conclude in Section PVTl with 
discussion on the application of these results to more gen- 
eral late-merger-ringdown models, such as the implicit 
rotating source model of Refs. [25|, [26[. We present a de- 
tailed description of our selection of equivalence classes 
of binaries in [Bj 



II. BUMPS IN NUMERICAL (3,2) MODES 

The first gravitational waveforms extracted from nu- 
merical simulations were the dominant (2, ±2) modes, 
whose early-inspiral behavior was expected to match the 
quadrupole radiation predicted by quasi-Newtonian and 
post-Newtonian theory. Once these had been shown to 
be robust and universal across codes [3, [|| , some groups 
turned their attention to the subdominant modes. An- 
alyzing the subdominant modes of equal-mass binaries, 
Buonanno et al. Q reported that an accurate fit of the 
(3, 2) mode for the ringdown stage of effective-one-body 
(EOB) waveforms requires the addition of the fundamen- 
tal (2, 2) quasinormal frequency. When Baker et al. [25[ 
looked at a set of mergers of nonspinning black-hole bi- 
naries with mass ratios in the range 1:1 to 6:1, they 
noted that one of the leading subdominant modes, (3, 2), 
showed an unusual bumpiness just after merger over a 
range of parameter space. This bumpiness manifested in 
both the frequency and amplitude, and appeared to per- 
sist with both increased resolution and extraction radius, 
thus constituting a robust pattern of excursions from the 
frequency separation dominating the i = m modes. More 
recent work by Kelly et al. [2(| shows the same anomaly 
in equal-mass binaries with non-preccssing spins (i.e., the 
spins are aligned/anti-aligned with the orbital angular 
momentum) . 

Examples of these more complicated waveform fea- 
tures are shown in Fig. [21 where we plot waveform fre- 
quency (top panel) and amplitude (bottom panel) of the 
measured (3, 2) mode for the merger of a nonspinning 
4:1 binary, as well as for the mergers of several other 
BH configurations with the same final dimcnsionless spin 
(ctf sa 0.475). We also mark the expected real QNM (2, 2) 
and (3, 2) frequencies, W22 and W32 for a Kerr black hole 
of this spin. From the time of peak amplitude (t = 
here) until the waveforms start to degrade around 60M 
later, the frequency seems to oscillate around one or other 
of these two QNM frequencies, rather than locking onto 
the higher W32, as for other modes. These oscillations ap- 
pear in the strain h and its time-derivatives; we choose to 
study strain-rate, h(t), waveforms, which we decompose 
into modes h^ tTn ){t) = \h(£ m )| exp(z(p(f m )), with instan- 
taneous frequencies (p(i, m )- 

We can model the more complicated ringdown wave- 
form features by expressing the (2, 2) mode as a pure 
QNM ringdown, and the measured (3, 2) mode as a lin- 
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FIG. 2: Post-merger frequency [top panel] and amplitude 
[bottom panel] of the numerically measured (3, 2) mode for 
a set of "4:l-equivalent" evolutions, resulting in a final black 
hole with dimcnsionless spin at « 0.475, matching that of 
a 4:1 nonspinning binary merger. The data sets have been 
shifted in time so that t = corresponds to peak amplitude 
of the dominant (2, 2) mode. The two dashed (black) hor- 
izontal lines in the top panel mark the fundamental QNM 
frequencies U22 (lower) and W32 (higher) for a Kerr hole of 
the same final spin. 



ear combination of QNM ringdowns: 

A 22 e i ^ 22t+522) , (1) 



^ model 
(2,2) 



t model 
n (3,2) 



Here o~£ m = ui£ m + i/r£ m is the full complex QNM fre- 
quency, and P32 = po exp(i£) is a constant complex- 
valued parameter indicating the mixing of the (2, 2) QNM 
mode into the measured (3, 2) mode. The modeled (3, 2) 
mode frequency and amplitude are then: 
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F(t) 
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where F(t) = 1 + 2e(t) cos(A R t + 5) + e(t) 2 , A R = oj 32 - 
W22, A/ = l/r 32 - l/r 22 , e(t) = p A 22 /A 32 exp(Ajt) = 
eo exp(Aji), and <5 = 5 32 - S22 — (. 

For a given mass and spin, the QNM frequencies, uj 22 
and W3 2 , and the damping times, t 22 and T32, are values 
known from black hole perturbation theory. Typically, 
t 22 Ri T32, so that A/ is somewhat smaller than Ah, 
allowing a beat-like effect to persist over several cycles. 
Fixing these leaves just two free parameters for the fre- 
quency: Eq = poA 22 / 'A32, the initial ratio of contributing 
amplitudes, and 5, the initial phase difference, as well as 
one more amplitude parameter, ^3 2 . 

Evidently, the characteristic shape of the modeled (3,2) 
mode frequency plots will depend on the relative magni- 
tude of the modal contributions: for eq C 1, the fre- 
quency will oscillate (approximately) sinusoidally about 
w 32 ; for Eq 1, the oscillation will be about w 22 ; for 
intermediate values, the oscillatory shape will be more 
complex. In the left panel of Fig. wc demonstrate 
these shapes for a Kerr hole of spin a>[ = 0.475, the same 
4:1 end-state spin as in Fig. [2] Similarly, the right panel 
shows the corresponding modal amplitude shape for the 
same end-state hole. Again, the most extreme bumps 
in amplitude occur when the (2, 2) and (3, 2) modes have 
comparable amplitude contributions (eo ~ 1). These the- 
oretical curves should be compared with the numerically 
measured mixing in Fig. [5J 



III. POSSIBLE CAUSES OF MODE-MIXING 

The bumpy features seen in the measured (3, 2) mode 
are a clear exception to the general rule that each angu- 
lar mode encodes sinusoidal waves with just one slowly 
evolving frequency component, the phenomenon we re- 
fer to as frequency separation. In Fig. [3J we showed 
that a combination of the fundamental (3, 2) and (2, 2) 
quasinormal-mode frequencies produces similar features. 
More generally, there are indications that such mixing 
occurs among other modes, especially other higher-order 
m = 2 modes, which likewise seem prone to coupling to 
the dominant mode. Here we ask the basic question: is 
this mode-mixing a fundamental property of the radia- 
tion, or some kind of an artifact, and if so, what kind? 

We consider various hypotheses to explain this mode- 
mixing effect violating our empirical frequency separation 
rule. The first, which we label physical mixing, is simply 
that the frequency separation rule does not physically 
hold to sufficiently high precision; that is we are are per- 
haps seeing a nonlinear effect in the radiation-generation 
process underlying the (3, 2) mode. Under this assump- 
tion, no choice of fixed or slowly evolving angular basis 
could be expected to yield the kind of frequency sep- 
aration we see in other cases. Near the merger where 
nonlinear physics is dominant, it is difficult to make any 
strong argument for expecting frequency separation. In- 
deed, we would be surprised to not find violations of this 
assumption as wc probe beyond the first few orders of 
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FIG. 3: Possible shapes from equations ((3j)-((4j) for the mod- 
eled (3, 2) mode frequency [top panel] and amplitude [bottom 
panel] resulting from a merger with an ai — 0.475 endpoint. 
All curves assume zero phase difference S, and overall am- 
plitude is arbitrarily scaled to unity at t = 0. As with fig- 
ure [2] the two dashed (black) horizontal lines in the top panel 
mark the fundamental QNM frequencies CJ22 (lower) and W32 
(higher) for a Kerr hole of the same final spin. 



magnitude in waveform precision. 

In the linear ringdown dynamics where this investiga- 
tion is focused, some degree of physical frequency sepa- 
ration can be expected, based on the separability of the 
Teukolsky equation, which describes small distortions of 
a stationary black-hole spacetime. The scale of physical 
linear mode-mixing can be quantified by careful consid- 
eration of quasinormal modes. 

The alternative hypothesis is that the mixing is an ar- 
tifact of our analysis, arising from choices that we make 
in setting up the angular mode decomposition. Perhaps 
our basis is not quite optimal, but we can find some other 
basis in which we more precisely recover frequency sep- 
aration. Indeed, given the freedom available in selecting 
such a representation, we have little grounds for suppos- 
ing that our first guess would be optimal. Here we con- 
sider two classes of choices in how to represent the space 
of gravitational radiation waveforms, which, in the full 



sense, has angular and retarded-time dimensions. 

The first choice we make is in how we define the 
spheres on which angular harmonic decomposition will 
be conducted. Within the structure of asymptotically 
flat spacetimes, gauge freedom in the choice of constant- 
rctarded-time spheres can yield a frequency-dependent 
mode-mixing effect in the decomposed waveforms. This 
ambiguity arises from the freedom to re-parameterize the 
proper-time coordinate, the so-called "supertranslations" 
subgroup of the Bondi-Metzner-Sachs gauge group for 
outgoing radiation. We describe this possibility of super- 
translation gauge mixing in more detail below. We may 
generally expect that mode-mixing of this sort will be 
most evident in the late merger, where wavelengths are 
shortest. 

The next choice we make is in choosing the family of 
angular basis functions on the extraction spheres. In this 
case, the mixing arises if our chosen family of modal ba- 
sis functions used for radiation extraction differs from 
the optimal one in which frequency separation is best 
approximated. It is common to apply a spin-weighted 
spherical-harmonic basis, but a different choice may be 
motivated for the ringdown signals. Indeed, the separa- 
tion of the Teukolsky equation is not achieved in a spin- 
weighted spherical harmonic basis, but in a spin- weighted 
sp/iermGW-harmonic basis. It has been suggested 0, [H[ 
that this difference explains the sort of waveform phe- 
nomena we consider, though this has not been demon- 
strated. We label this effect angular-basis mixing. 

In the next subsections, we consider these possible mix- 
ing effects in detail, preparing for a quantitative study of 
the evidence for these effects in numerical data in Sec- 
tion IVl 



A. Gauge Effects 

To understand the effect we are calling supertransla- 
tion gauge mixing, we must make a brief detour to de- 
scribe the gauge freedom in the representation of an out- 
going radiation field approaching future null infinity in an 
asymptotically flat spacetime. Consider such a spacetime 
in standard retarded-time coordinates {it, r, 9, 0}. Scaled 
by r, the outgoing radiation field propagates outward on 
null rays labeled by u, 9, and (f>. Each polarization com- 
ponent can thus be described by a function of these vari- 
ables. The Bondi-Metzner-Sachs (BMS) [H SI group 
describes gauge transformations among these variables 
of the form 



9' = 0'(M), 4>' = </>'( 



u' = K{9, <P) (u - a(9, </>)) 



where (#,</>) — > (9',</)') is a conformal transformation on 
a constant-u sphere with conformal factor K . 

For concreteness in the context of numerical relativity 
simulations, we note that it is common to make these 
gauge choices by specifying an "extraction sphere" lo- 
cated sufficiently far from the source where radiation field 
calculations arc realized. The effect of one class of BMS 
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transformations, amounting to rotations of the extraction 
sphere, has been identified as an important concern when 
the choice of axis is not fixed by symmetry [3414391 ] . How- 
ever, the simulations in this study involve nonprecessing 
mergers, with no ambiguity in defining the orientation of 
the extraction sphere. 

But what happens if we make a small radial perturba- 
tion of the extraction sphere? It is clear that sufficiently 
small distortions of larger extraction spheres would have 
negligible impact on the intrinsic geometry of the sphere. 
The gauge effects of such distortions are described by a 
subset of the BMS transformations, known as supertrans- 
lations, with 8' = 6,<j>' = <f>, and K = 1. 

Now consider the effect of a supertranslation on a grav- 
itational waveform i/)(u, 9, </>). Here we will make the ad- 
ditional assumption that a(0, 4>) is sufficiently small that 
we can approximate the effect of the supertranslation by 



a{e,cj>)-^{u,6,4>), (5) 



and we can expand the supertranslation in terms of 
(scalar) spherical harmonics: 



(6) 



LM 



Then from ([5]), the measured radiation modes will be 
perturbed as follows: 

ll>e.m{v!) ~ 1plm{u) + ^ Ci m t>m'-Q-i>t'm'(v>), (7) 



where 



E 6 

LM 

E 5 



r L -21 e> 



LM 



LM 



(2L+l)(2f + 1) 
4tt(2£+ 1) 



1/2 



x (L, 0, £', 2\£, 2) (L, M, £', m'\£, m) . (8) 

In this paper we focus on mixing from the dominant 
mode, (£' = 2, m! = 2), with another m = 2 mode, fixing 
these values. For these cases the Clebsch-Gordan selec- 
tion rules require that M — and £ — 2 < L < £ + 2. 
Then our mixing coefficient takes the form 

1/2 



C^222 — ^2 ^LO 



5 (2L + 1) 
47(27+1) 



(i,0,2,2|^,2) z . (9) 



For example, complete expansions for £ = 3 and £ = 4 
would yield 
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FIG. 4: Real and imaginary parts of p ga ug C ,32 and p gaU ge,42 
(for &20 = 1) for post-merger Kerr BHs of dimensionless spin 
«f (and corresponding fundamental (2, 2) QNM frequency 

MfO-22). 



The shape of the distorted extraction sphere is deter- 
mined by the coefficients : for real a, we need the 
&lo also t° be real. The reality of the Clebsch-Gordan 
coefficients then implies that CV222 is also real. 

The other ingredient in the waveform-mode perturba- 
tion (JT)) is the derivative w.r.t. u on the right-hand side: 



d u ^ e , m ,{u)=d u (A{u)e 1 ^ 




'(«)• 



After merger, the effective coefficient \ A/A + iipj will 
asymptote to a constant complex number: 




1 



+ iuii' m > = ic 



This implies a simple, QNM-driven leakage from the 
(2, 2) mode into higher-^ modes. Collecting terms, and 
working with the strain-rate h, during ringdown we have 



J, gauge , 
"■(£,2) ~ ft (^ 2 



iCl222 &22 h(2,2) 
h(l,2) + Pgauge,^2 ^(2,2) ■ 



(10) 



In Fig. [4J we show the real and complex parts of the 
leakage parameters p ga uge,32 and p ga uge,42 for the sweep 
of end-state spins at, assuming an unchanging scaling 
620 = 1 (and all other b^o = 0). The value of 6 20 is 
not physical, but gauge, and may differ between any two 
waveform determinations. The most important property 
we note is that the BMS leakage coefficients are nearly 
pure-imaginary at any fixed 620 and any spin af. 
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B. Angular Basis Effects 

Another possible path to mixing arises from consid- 
ering what quasinormal-mode (QNM) frequencies actu- 
ally represent. QNMs were originally discovered in nu- 
merical black hole scattering studies [4(| El[ and even- 
tually understood as a key feature of the perturbation 
theory of Kerr black holes [42| ■ In developing this theory, 
Tcukolsky worked with a background Kerr black hole in a 
very specific coordinate system due to Boyer & Lindquist 

A perturbed Kerr black hole will ring down to qui- 
escence through the emission of gravitational waves. 
These waves will have characteristic frequencies u>£ m 
and damping times r^ m given by the hole's QNM spec- 
trum. 2 While the primary aim of QNM analysis is 
to determine the set of allowed complex frequencies 
0( m = u!i rn + i/Tt m , these frequencies are tied to the 
radial and angular eigenfunctions arising from the sep- 
aration of the perturbation equations. These angular 
eigenfunctions are the spin-weighted spheroidal harmon- 
ics, _ 2 ;yr(Afacr; #,</>) = - 2 S]T {Maa; cos 6)e im < t '. 3 

Numerical waveform extraction from binary mergers, 
on the other hand, typically decomposes the waveforms 
onto the more generally motivated basis of spin- weighted 
spherical harmonics ^ 2 Y^ n {6 1 </>), which correspond to a 
spheroidal harmonic basis with Maa = 0: - 2 Y e m (6, (f>) = 
-2^ m (°; 0. 0) S3- Buonanno et al. Q demonstrated that 
using a spherical harmonic basis will necessarily result in 
mixing of (£,m) and (£',m) quasinormal modes. With- 
out an obvious nontrivial choice for Maa that applies 
at all times, for all modes, over the course of the evolv- 
ing simulation, decomposing with Maa — > seems a 
natural choice. Here we consider an alternative choice, 
Maa — > MfQ!fcr22, hoping to limit much of the mode 
mixing. Using this basis requires knowing the final Kerr 
state (Mf , at) of the merger before the decomposition can 
be applied, and the additional task of numerically com- 
puting the basis functions (see Appendix |A"|) . Still this 
basis is not optimal for the subdominant modes. This un- 
avoidable sub-optimality is discussed further in the next 
subsection. The distinction between the spheroidal and 
spherical harmonics may be expected to yield the appear- 
ance of mode-mixing in the numerical waveform results 
even if we have eliminated the gauge freedom noted in 
the last section by optimal correspondence with a suit- 
ably perturbed Boyer-Lindquist coordinate system. 

To estimate the apparent mode-mixing from this ba- 
sis mismatch, we can calculate the overlaps between the 
spheroidal harmonics (for a particular Maa) and the 



1 Teukolsky theory can be reformulated on other backgrounds; see, 
e.g., 0. 

2 We omit the principal quantum number n, assuming that we are 
dealing with the slowest-damped fundamental (n = 0) QNM. 

3 Here we use the symbol a to denote a generic complex frequency. 
(T£ m is a specific eigenvalue of the Kerr background. 



spherical harmonics. That is, we want to know the coef- 
ficients S(ii m in 

oo 

_ 2 yF(Maa;e,<f>) = £ s t , lm _ 2 Y^[B^). (11) 

i'=i 

We describe our calculation of the - 2 y™ in Ap- 
pcndix[XJ To determine the overlaps s^ m , we decompose 
the properly normalized spheroidal harmonic against the 
spherical harmonics in the usual way: 

sti m =jdn - 2 yr(Mia { a 22 ;e, <j>) -2^(0, 0)* 

= J dx- 2 Sr(M f a f a 2 2;x)- 2 SP{0;x)*. 

Now consider the idealized case where a physical ring- 
down signal is the simple combination of the fundamen- 
tal (2,2), (3,2), and (4,2) quasinormal modes (we omit 
Maa arguments for brevity): 

4 

h(t, r, e, <!>) = Hn(t, r) - 2 yl{e,<t>) 

i 

4 

» J>* a (*»»0_ 2 i2(M). (12) 

l> 

If we make the reasonable assumption that mixing I 7^ I' 
products can be ignored for subdominant modes, then 
the measured spherical harmonic ringdown modes are ap- 
proximately: 

tf2%)(t,r) « s 2 , 22 H 22 (t,r), 

tfe%&r) ~ '*aHa(t,r) + p hasis ,i 2 h 2 , 2 {t,r). (13) 
Here, the mixing coefficients arc 

— Si " 22 flA\ 
Pbasis,£2 = ■ U^j 

S2'22 

In Fig. [5J we plot the coefficients pbasis,32 and /9basis,42, 
evaluated at Maa = M{ata 22 , where a 22 is the funda- 
mental QNM frequency of the (2, 2) mode for a Kerr hole 
of mass Mf and dimcnsionlcss spin af. Note that (a) 
there is no ambiguity in overall scale for these coefficients 
(unlike the BMS-derived coefficients of the last section) , 
and (b) they are strongly real-dominated. 



C. Physical mixing 

The discussion above exposes artifacts that arise from 
waveform decomposition using ordinary spin-weighted 
spherical harmonic functions. Here we ask whether, even 
with extraction spheres in the Boyer-Lindquist gauge, an- 
other decomposition using spin- weighted spheroidal har- 
monic functions can avoid mode-mixing. 
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FIG. 5: Real and imaginary parts of the mixing coefficients 
pbaaia,32 and pbasis,42 for post-merger Kerr BHs of dimension- 
less spin Qf (and corresponding fundamental (2, 2) QNM fre- 
quency Aff<722). At the zero-spin limit at — ¥ 0, the leakage 
vanishes. 



The question is non-trivial. Although each leading- 
order quasinormal ringdown mode exhibits angular de- 
pendence described by some kind of spin- weighted spher- 
ical harmonic angular function, they are not mutually 
given by the same kind of spin-weighted spherical har- 
monic angular functions, since each has its own distinct 
quasinormal frequency agm, and consequently a distinct 
preferred basis as labeled by Maa = M[a[<Tg m . We must 
choose some particular orthonormal basis for the decom- 
position, and that basis cannot be simultaneously opti- 
mal for each mode. 

That the spheroidal harmonics associated with differ- 
ent QNM frequencies are not perfectly orthogonal has 
been demonstrated for high-spin Kerr holes by Berti et 
al. [45i | . To quantify this for a general end-state spin at, 
we define new overlaps tin, between spheroidal harmon- 
ics associated with different m = 2 QNM frequencies: 



tvt = jdVl _23?(Mf«f<T«) - 2 y t ? (MfOfO-^a)*. (15) 

Figure [6] shows the magnitude of these overlaps for I = 2 
(upper panel) and I = 3 (lower panel), and several val- 
ues of £'. From these plots, we see that the spheroidal 
harmonics for different M aa are not orthogonal, but 
show mixing by as much as 4% for high spins (though 
the maximum overlaps occur at sub-maximal spins, as 
noted by [45(). The overlaps are also greatest for "near- 
est neighbor" modes: I = £' ± 1. For example, if we 
decomposed a waveform, including a non-trivial (2,2) 
QNM, in the spheroidal basis corresponding to the (3, 2) 
mode ringdown frequency, then the corresponding curve 
in Fig. [5] would represent a mixing coefficient analogous 




FIG. 6: [Top panel] Magnitude of overlap t(2 between -2^21 
evaluated at Maa — M{Ctf(T22 and ~2y'i, evaluated at 
Maa — MfQf<Tf2. [Bottom panel] Same for overlap tt3 be- 
tween -2^3 and -iy\. By definition, |t32| = [^23 1 - 



to those in the previous subsections. There is no choice of 
orthonormal basis that will avoid all such mode mixing. 
In this sense, the angular non-orthogonality of the quasi- 
normal mode implies a form of physical mode-mixing, 
meaning that we can not perfectly isolate the QNM fre- 
quencies by any choice of angular basis. 

Fortunately it seems that the most evident mixing in- 
volves the dominant (2,2) mode frequency bleeding into 
higher-i! modes. With that assumption we may still elim- 
inate most physical mixing by choosing the basis compat- 
ible with this dominant quasinormal mode. If we decom- 
pose with the basis labeled by Maa = MfafO-22 then 
the orthogonality of this particular basis will completely 
prevent the (2, 2) quasinormal mode from mixing in to 
any other decomposed modal waveform component. In 
this way we can eliminate any "physical mixing" of the 
particular form described in Section [TTJ Mixing among 
subdominant modes, or mixing of subdominant modes 
into the decomposed (2, 2) waveform component will still 
occur at some level, but this is a smaller effect, which we 
do not focus on in this paper. 
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IV. SIMULATIONS 



B. Waveform Extraction 



To investigate the mixing in a systematic way, we have 
surveyed several existing simulations of aligned-spin bi- 
naries, as well as carrying out new short simulations 
with the Goddard Hahndol evolution code. We choose 
our new black-hole binary (BHB) configurations in sev- 
eral groups of "merger-equivalent" classes, as described 
in [B] The initial parameters for all these simulations, 
old and new, are presented in Table |U In Fig. [JJ we 
show the distribution of these configurations as plots in 
the two-dimensional configuration-spaces {01,02} and 
{ai,g}, where q = Mi/M% > 1 is the mass ratio, and 
a a = Sa I is the dimcnsionlcss spin parameter of hole 
A, with physical values restricted to a a € [—1,1]. Many 
of the longer and higher-resolution evolutions have ap- 
peared in previous publications [25|, [26| . Since our pri- 
mary interest here is strictly in the late-merger regime, 
newer evolutions begin only a few orbits before merger. 



Numerics 



The initial momenta of the newer evolutions were cho- 
sen by integrating the post-Newtonian equations of mo- 
tion, as outlined in [3a . |47| , with spin contributions to 
the Hamiltonian adapted from [48145 2j , and the flux from 
[53|. Note that we did not attempt to reduce the eccen- 
tricity through tuning the initial momenta. 

The new evolutions use the Hahndol codepaired with 
the "Curie" release of the Einstein Toolkit 1541 . incorpo- 
rating the Cactus Computational Toolkit [55[ and the 
Carpet mesh- refinement driver (56j . 

In all cases, the initial data are of the standard Brandt- 
Briigmann puncture type 57 1, using the Bowen-York (58j 
prescription for extrinsic curvature that exactly satis- 
fies the momentum constraint. We solve the remaining 
Hamiltonian constraint using the TwoPunctures spectral 
code (H. 



To evolve these initial data, we employ the BSS- 
NOK 3+1 decomposition of Einstein's vacuum equa- 
tions [60| - [62j . with the alternative conformal variable 
W = e~ 2< ^ suggested in fo3l - l65j . constraint-damping 
terms suggested in [66| . and the dissipation terms sug- 
gested in [67], H|[ . Our gauge conditions are the specific 
1+log lapse and Gamma-driver shift described in [o9| . 
which constitute a variant of the now-standard "moving 
punctures" approach 0, Q . Our spatial derivatives use 
sixth-order-accurate differencing stencils, with the excep- 
tion of advection derivatives, where we use fifth-order- 
accurate mesh-adapted differencing (MAD) [7(J. Our 
time- integration is performed with a fourth-order Rungc- 
Kutta algorithm. 



We extract the gravitational waveforms from the sim- 
ulations through the radiative Weyl scalar -04 [28| . This 
is evaluated throughout the grid, and interpolated onto 
a set of coordinate spheres at extraction radii r € 
[40M, 90M]. Over each sphere, the interpolant is inte- 
grated against the set of spin- weighted spherical harmon- 
ics _ 2 *7™(M), up to I = 5. 

In the extraction region, the grid spacing is between 
M/2 and 2M, depending on the central resolution of the 
simulation. This is generally too coarse to resolve higher- 
frequency (and higher-m) modes with accuracy. Even for 
the dominant, relatively low-frequency, (2, ±2) modes, 
dissipation effects are visible that spoil the 1/r extrap- 
olation near and after merger. For this reason, we have 
used an r-extrapolation scheme that includes an explicit 
dissipative term in the amplitude of each mode: 

At m (r) = a_ir + a + a 2 r~ 2 , ip im (r) = b + b 2 r~ 2 . 

We have found this extrapolation procedure to be robust 
only for the higher-resolution simulations in this paper. 

As a result, a waveform-derived quantity / will have 
errors due to finite extraction radius and finite resolu- 
tion. For this paper, we make a very conservative error 
estimate by adding uncertainties linearly: 

A/ = A,,./ + A h f. 

For the finite-r error, we assume an uncertainty equal 
to the difference between the coefficient from the r- 
cxtrapolated highest-resolution data and that measured 
from the largest finite-r data at the same resolution. For 
finite-resolution error, we use the difference between the 
same-extraction-radius data at the coarse and fine reso- 
lutions as our estimate of the error in the fine-resolution 
result. For many configurations, we only have a single 
resolution available and the r-extrapolation is not reli- 
able at this resolution. For these, we adopt a conser- 
vative overall error estimate by taking the average error 
from comparable two-resolution configurations 4 and mul- 
tiplied it by 1.5. For amplitudes, this was a relative error, 
while for phase measurements, it was the absolute error. 



V. ANALYSIS OF WAVEFORMS 

Using the ringdown data from all the simulations in 
Section IIV1 we performed least-squares fits to the real 
part of the strain-rate (2, 2) and (3, 2) waveforms, using 
the forms of equations ([I])-©- Our fit is over the win- 
dow t G [20,55], where t = is the time of peak (2,2) 



4 By "comparable" , we mean configurations that used the same 
numerical executable and grid structure, and whose lower- 
resolution version matched that of the single-resolution configu- 
rations. 
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TABLE I: Physical and numerical parameters of the initial data for all the runs presented. mi lP and 7712, p are the bare puncture 
masses of the two pre-merger holes, ro is the initial coordinate separation, while Pot and Po r are the initial transverse and 
radial components of the Bowen-York linear momentum. Madm is the total energy of the initial data, while the total infinite- 
separation mass of the system is estimated by the sum of the initial ADM masses of the individual holes We have found 
that for all cases here, this differs from the sum of apparent-horizon masses (calculated at times between t = 100 and 200), by 
less than a tenth of a percent. 



run name 




ITl2,f 


Si. 


S2z 


ro 


Pot 
(xlO 2 ) 


Por 

(xlO 4 ) 


Madm 


5Zi ^ADM,i 


X1JJU 


0.301805 


0.301805 


0.2 


0.2 


8.20 


10.32 


0.00 


0.988459 


1.000804 


Xl.uu 


0.454575 


0.454575 


0.1 


0.1 


10.21 


9.25 


9.17 


0.99223 


1.002768 


X1.00 


0.487231 


0.487231 


0.0 


0.0 


11.00 


9.01 


7.09 


0.990514 


1.000050 


X1JJD 


0.301805 


0.301805 


0.2 


-0.2 


11.00 


9.01 


7.09 


0.990024 


0.999222 


XI . 5_00 


0.581359 


0.380645 


0.0 


0.0 


7.12 


11.75 


29.17 


0.987252 


1.000000 


XI. 75.00 


0.619237 


0.345598 


0.0 


0.0 


7.42 


11.01 


24.10 


0.988129 


1.000000 


X2.00 


0.649344 


0.314904 


0.0 


0.0 


7.00 


11.00 


0.00 


0.987939 


1.000000 


X2.DU 


0.648662 


0.265507 


-0.066666667 


0.066666667 


10.00 


8.52 


7.63 


0.990951 


1.000009 


X2.5.00 


0.699349 


0.269501 


0.0 


0.0 


7.40 


9.79 


20.53 


0.989664 


1.000000 


X3.00 


0.738687 


0.237505 


0.0 


0.0 


8.88 


7.88 


8.96 


0.991673 


1.000000 


X4.00 


0.790000 


0.189000 


0.0 


0.0 


8.47 


6.96 


0.00 


0.992912 


1.000310 


X5JJ0 


0.822007 


0.157080 


0.065083333 


0.0 


8.68 


5.91 


4.88 


0.993733 


1.000000 


X3_d0 


0.731667 


0.237705 


-0.087566063 


0.0 


9.06 


7.84 


8.76 


0.99187 


1.000000 


X2_d0 


0.587677 


0.317821 


-0.210380889 


0.0 


8.44 


9.93 


16.53 


0.989967 


1.000000 


X1.DD 


0.390411 


0.390411 


-0.159125 


-0.159125 


11.98 


8.84 


1.20 


0.990453 


0.998786 


X5.00 


0.824897 


0.157031 


0.0 


0.0 


8.67 


5.97 


5.85 


0.993827 


1.000000 


X6.00 


0.848615 


0.133064 


0.0 


0.0 


7.55 


5.84 


6.94219 


0.994008 


1.000000 


X5.D0 


0.822405 


0.156318 


-0.052232639 


0.0 


8.09 


6.32 


7.00085 


0.993556 


1.000000 


X4.D0 


0.778549 


0.188766 


-0.1213184 


0.0 


8.57 


7.04 


7.78076 


0.992926 


1.000000 


X3.D0 


0.692530 


0.237756 


-0.21614625 


0.0 


9.17 


7.93 


8.80172 


0.99219 


1.000000 


X2.DD 


0.531347 


0.260245 


-0.277766667 


-0.069441667 


10.72 


8.56 


7.81844 


0.992008 


1.000000 
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FIG. 7: Simulated configurations from Table[T] represented as points in two-dimensional {qi, Q2}-space [left panel] and {a lt q}- 
space [right panel]. Note that some points are associated with multiple simulations. 
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mode amplitude. By starting 20M after peak amplitude, 
we ensure that we are in the linear ringdown regime; 
by stopping at 55M, we avoid the low-amplitude degra- 
dation seen in late-ringdown waveforms. As the tabu- 
lated version of the results would be excessively long, we 
present our raw results purely graphically. 

We begin by showing the nature of the complex numer- 
ical "leakage parameter" derived from the ratio of fitted 
parameters from the measured (2, 2) and (3, 2) modes 
during ringdown, using ([T])-([2]): 



Pnum,32 



p 32 A 22 e^ t+s ^ 



(16) 



Figure [5] shows the real and imaginary parts of this leak- 
age for all configurations presented in this paper, as a 
function of the dimcnsionlcss spin ai of the post-merger 
hole. 



A. Comparing Hypotheses 

In Section Mil we discussed two possible causes for 
mode-mixing effects of the form 



\ model 
n (i,2) 



A l2 e 



pe 2 A 22 e 



i{<722t + &22) 



(17) 



described in Section [II] If the mixing is caused by BMS 
supcrtranslation gauge ambiguity, then we would expect 
nearly pure imaginary pg 2 . On the other hand, if the 
mixing derives from the distinction between spheroidal 
and spherical harmonic angular functions, then we expect 
predominantly real pg 2 of a quantified size. In Fig. [5] we 
see that the argument of p num ^ 2 is close to zero, within 
error bars for most cases, making p n um,32 predominantly 
real, consistent with the spheroidal harmonic hypothesis. 
The largest deviations from zero are also those with the 
largest uncertainties arising from the QNM fit process. 

The analysis of Section [ill AI suggests that a change of 
supcrtranslation gauge would give rise to mode-mixing 
coefficients with a numerically significant imaginary part 
in the measured waveform. Since the imaginary part of 
Pnum,32 is so small, any supertranslation gauge effects are 
negligible at the level of interest here. We can estimate 
the degree of gauge constraint implied by a null mea- 
surement of this effect. From the bottom panel of Fig. [51 
one sees that the imaginary component of the mixing co- 
efficient /9 n um,32 is constrained to values within ±0.02 in 
almost all cases. If we generously assumed that all of this 
imaginary mixing was caused by gauge distortion of the 
extraction sphere, by comparison with Fig. [U we would 
conclude that the amplitude of the distortion (620 specif- 
ically) would be have to be smaller than about 0.2M, 
suggesting a remarkable level of supertranslation gauge 
optimality in these simulations. 
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FIG. 8: Real [top panel] and imaginary [bottom panel] parts 
of complex leakage parameter pnum.32 (|16p . shown as a func- 
tion of at , the dimensionless spin of the post-merger hole for 
all runs, indicating that in almost all cases, the leakage pa- 
rameter is predominantly real. 



B. Testing the Spheroidal Leakage Model 

We have seen that the numerical results for the com- 
plex argument of p 32 are consistent with the spherical- 
spheroidal mixing hypothesis, but this hypothesis also 
makes quantitative predictions for the magnitude | P32 1 - 
In the top panel of Fig. [SJ we plot the magnitude of 
Pnum,32 as a function of af. We overlay these points 
with the magnitude of the leakage coefficients pbasis,32 
(Q3J) plotted in Fig. [5] (blue solid curve). From the close 
fit, it appears that the leakage is in fact dominated by 
this spheroidal/spherical harmonic mismatch. That is, 
even though the post-merger background coordinate sys- 
tem should not be expected to closely resemble the Kerr- 
Boyer-Lindquist slicing assumed by Teukolsky's pertur- 
bative work, nevertheless, this expected warping is not 
as important as our choice of harmonic basis functions. 

The bottom panel of Fig. [H] shows the complex ampli- 
tude of the equivalent parameter p n um,42 governing the 
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FIG. 9: [Top panel] Amplitude of measured (2, 2) -»■ (3, 2) 
leakage parameter p numj 32 as a function of at, the dimension- 
less spin of the post-merger hole for all runs. The curve is the 
theoretical ratio pbasis,32 due to spheroidal-spherical leakage. 
[Bottom panel] Same for (2, 2) — y (4, 2) leakage ratio p n um,42. 



leakage of the (2, 2) mode into the measured (4, 2) mode. 
Although this is also consistent with expectations from 
angular-basis mixing (blue solid curve), the relative er- 
rors swamp the numerical data, and higher-resolution nu- 
merics will be needed to establish the relation unambigu- 
ously. 



C. Finding the Residual (3, 2) Mode Amplitude 

If we regard the measured (3, 2) mode as the combi- 
nation of a "true" (3,2) mode A32 exp 1(032^ + fe) and 
a piece of the (2,2) mode, we may ask whether we can 
model the residual (3, 2) contribution. When looking at 
the entire suite of simulations, it is difficult to see a dis- 
tinct pattern in these true (3, 2) amplitudes. However, it 
is instructive to carry out a particular slice in configura- 
tion space. 

In Fig. [TU1 we show a subset of the (3, 2) amplitudes 



0.02 



0.015 




0.005 



FIG. 10: Corrected amplitude A32 of the (3, 2) strain-rate 
mode for nonspinning binaries with mass ratio q € {1.0, 6.0}. 
The solid (green) curve is a fit to these points of the functional 
form (fig)) . 



formed by the mergers of nonspinning binaries, with 
mass ratio q = M1/M2 E {1.0,6.0}. Error bars in 
this plot have been estimated in the same way as for 
Fig. Clearly the high-g behavior seems to decay to 
some constant amplitude, while there is some local min- 
imum around rj = 0.21 (between q = 2 and q = 2.25), in- 
dicating that perhaps at this mass-ratio, the (3, 2) QNM 
is hardly excited at all. 

We also present an empirical fit to this data of the 
functional form 



A 32 (v) = \f {a - be-V«) 2 + c 2 



(18) 



where the parameters take the values a = 0.0147±0.0002, 
b = 1.5 ± 0.1, c = 0.0026 ± 0.0002, and A = 0.98 ± 0.02. 
Fits of this form are expected to be useful for generating 
merger template waveforms for the subdominant modes. 



VI. DISCUSSION 

In this paper, we have investigated "bumps" measured 
in the merger-ringdown portion of certain gravitational- 
radiation angular waveform modes from the numeri- 
cal simulation of the coalescence of black-hole binaries 
(BHBs). These bumpy modes appear to contain signifi- 
cant contributions from the dominant (2, 2) mode, indi- 
cating some kind of mode-mixing at work. 

We have considered three classes of effects that may 
contribute to mode-mixing in numerically extracted and 
decomposed merger-ringdown waveforms. These are: 
gauge effects, arising from supertranslation gauge free- 
dom for outgoing radiation in general asymptotically 
flat spacetimes (see Sec. IIII A"|) : angular-basis effects, re- 
lating to a choice between spin-weighted spherical or 
quasinormal-modc-adapted spheroidal harmonic bases 
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(Sec. IIII B|) ; and physical quasinormal-modc mixing ef- 
fects that are independent of any representation changes 
(Sec, lira . 

We have identified and analyzed the measured mode- 
mixing bumps in the most prominent of the bumpy grav- 
itational waveform modes modes — £ = 3,m = 2 — mea- 
sured from a set of numerical evolutions of aligned-spin 
BHB mergers. Our analysis has allowed us to distinguish 
between the contributions of our three mode-mixing ef- 
fects. We find that the angular-basis effects dominate. 
Although other kinds of effects may be present - like the 
frequency-dependent gauge supertranslations discussed 
in Section fill Al - they cannot be seen clearly here with 
the level of accuracy available from our current simula- 
tions. 

In this way our analysis further codifies the results 
from the ringdown stage of the aligned-spin mergers. 
This was originally prompted by our work on a multi- 
mode waveform model based on the implicit rotating 
source (IRS) picture of black- hole merger [25|,[26[. In this 
model, the dominant and leading sub-dominant wave- 
form modes from binary mergers were seen to share a 
common rotational phase, with a corresponding rota- 
tional frequency that increased monotonically through 
inspiral and merger, plateauing during ringdown. The 
corresponding mode amplitudes could be modeled by a 
simple, few-parameter functional form that depends on 
the frequency function, with a single well-defined peak. 
Attempting to extend this to the (3, 2) mode proved prob- 
lematic, as the measured mode was no longer monotonic 
in frequency, or single-peaked in amplitude. 

More broadly, we expect our results to provide guid- 
ance in the ongoing effort of combining results of ana- 
lytic and numerical relativity studies toward the goal of a 
fully developed family of efficient and accurate black- hole 
merger waveforms. Because the comparison of waveform 
models is typically conducted mode by mode in decom- 
posed form, the issues we have studied may lead to un- 
necessarily spurious features in particular waveform rep- 
resentations. 

We estimate, for instance, that supertranslation gauge 
changes that would effectively distort the shape of ar- 
bitrarily large waveform-extraction spheres on scales of 
order Mf or smaller would be sufficient to qualitatively 
influence the mode-mixing features focused on in this 
study. The absence of such effects is itself intriguing, 
suggesting that we have achieved nearly optimal choice 
of supertranslation gauge. Our near-optimal spheroidal 
harmonic basis is consistent with quasinormal-mode dis- 
tortions of Kerr space-time in the Boycr-Lindquist co- 
ordinate system. That we see negligible supertransla- 
tion mode-mixing suggests that the outer regions of our 
numerical space-times asymptotically approach distorted 
Kerr in Boyer-Lindquist coordinates faster (in powers 
of 1/r) than the asymptotic approach to perturbed 
Minkowski spacetimc. This seems plausible, based on 
our choice of numerical gauge, which approximates max- 
imal time slicing and V = spatial coordinates. The 



latter condition will yield spatially isotropic coordinates 
where possible. 

Nonetheless, it seems that we have been lucky to stum- 
ble onto a near-optimal representation as other incompat- 
ible gauge choices may also be reasonable in the numer- 
ical simulation context. In continued pursuit of higher- 
precision waveform comparisons and higher-fidelity ana- 
lytic models (see, e.g., the NR-AR project [7l|). we ex- 
pect such considerations to grow in significance. (They 
may also be crucial in studies of how the pre-merger 
BHB configuration is encoded in the relative amplitude 
of different quasinormal modes during ringdown; see, 
e -g- [z3-) Similarly we find that physical mode-mixing 
among the quasinormal modes will prevent any orthonor- 
mal representation from fully separating frequencies at 
sufficiently high precision. 

For simulations similar to ours, where gauge and phys- 
ical mixing effects remain small, and the primary source 
of mixing involves the (2,2) mode, our results suggest 
that decomposition with a spheroidal harmonic basis 
{ -2y£ l (MfCtf<j22', </>)} may be close to an optimal ba- 
sis for achieving modal frequency separation, and thus 
nearly beat-free waveforms. 
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Appendix A: Calculating Spheroidal Harmonics 

As there are no closed- form solutions for the -2^™, 
we must proceed numerically. While setting up the 
popular continued-fraction method for computing the 
quasinormal-mode (QNM) frequencies of a Kerr black 
hole, Leaver (73j presents the following power-series ex- 
pansion for the polar-angle function _2<S'|"'(M aa\ cos 6), 
due originally to Baber & Hasse [74| (we specialize here 
to s = -2): 

- 2 SP(Maa;x) =e Maax {l + x) |m+2|/2 (l - x) |m " 2|/2 

OO 

x^a„(l + x)", (Al) 

n=0 

where the expansion coefficients a n are determined up 
to an overall scaling — the value of ao — by the same 
recurrence relations that yield the QNM frequencies. For 
our desired Kerr spin af , we first determine the (complex) 
fundamental QNM frequency of the (2, 2) mode, Mf022- 
Next, assuming ao = 1, we use the recurrence relations 
from (731 ] to determine the a n (in practice, we truncate 
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the series at n = 14). Requiring that 

J dx\. 2 Sl n (M iai a 22 ;x)\ 2 = 1 

then fixes do, supplying the correct normalization of the 
a n . 

Appendix B: Kerr-Equivalent Black-Hole Binaries 

The end-point of any merger of BHBs in vacuum is 
expected to be a single Kerr black hole, parameterized by 
two numbers, the mass Mf and spin angular momentum 
Sf = otfMf . These should satisfy the global conservation 
rules: 



Mf — MadM — ^radi 
<Sf = ^ADM — ^rad7 



(Bl) 
(B2) 



where Madm and Jadm are the Arnowitt-Deser-Misner 
(ADM) energy and total angular momentum of the ini- 
tial data, and -E^ad and J ra <i are the energy and angular 
momentum emitted in gravitational radiation during the 
course of the evolution. 

Fixing the initial separation of the binary, and taking 
its total mass to be M = Mi + AI 2 (> Madm for any 
finite initial separation), and assuming zero eccentricity, 
the black-hole binary will have seven free parameters: 
{q, S\,S 2 }, where q = M\jM 2 > 1 is the mass ratio, 
and Sa are the spin angular momentum vectors of the 
two holes. However, the end-state has just two param- 
eters: {M{,Si}, so there must be a large degeneracy in 
the initial parameters. 

Viewing the BHB coalescence as a kind of simple par- 
ticle interaction, Boyle et al. [lH used symmetry argu- 
ments to restrict the possible end-states of the BHB 
merger. This is the basis of end-state formulae by Tichy 
& Marronetti [76|. Other models have been developed 
by Buonanno et al. Lousto et al. jz^], Barausse & 

Rezzolla [zD, H3| and others. 



In the case of initially orbit-aligned spins, the initial 
parameter space is three-dimensional: {q, S±, S 2 }. We 
use the simplest applicable formulas for the achieved end- 
state for an aligned-spin system. The end-state mass 
formula we take from Eq. (5) of [78j |: 

M f = 1 - rjEisco - E 2 if - E 3 if 

- 7— — 77 \Es{a 2 + q 2 ai) + Eg{l - q){a 2 - qa.\) 

(l + qy 

+E A (a 2 + q ai ) 2 +E D (a 2 -qa 1 ) 2 ], (B3) 

where r\ = MiM 2 /(Mi+M 2 ) 2 = q/(l+q) 2 is the symmet- 
ric mass ratio of the binary, and the fitting parameters 
are: 



v i ^-unifvwfM i (g(l + 2 g )q 1 + (2 + g ) Q2 ) 

ElSC ° =1 -~ + ' 103803?7 + 36V3(l + g)2 

5(gai - a 2 ) 2 



162V2(l + g) 2 ' 
E 2 =0.341, E 3 = 0.522, E s = 0.673, 
E s = - 0.36, E A = -0.014, E D = 0.26. 

For the final spin, one model with just enough com- 
plexity for our data sets here was given by [7^, |8(j 5 : 



,^2 



Qjf = a + siria^ + s 5 r] 2 a+t a r]a + 2V3r]+t 2 r] 2 +t 3 rj 3 , (B4) 

where a = (q 2 ai + a 2 )/(q 2 + 1) and the coefficients 
{s4,s 5l t 0l t 2l t 3 } are: 



s 4 = -0.1229 ± 0.0075, s 5 = 0.4537 ± 0.1463, 
t = -2.8904 ± 0.0359, t 2 = -3.5171 ± 0.1210, 
t 3 = 2.5763 ±0.4833. 

Using these formulae, we have constructed a set of con- 
figurations, which we present in Table HI grouped by final 
Kerr spin. 
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